inp<-
read.csv("primates-cr-registros.csv", sep = "\t")
inp %>%
dplyr::select(family,
species,stateProvince,
locality,eventDate) %>%
DT::datatable(
colnames = c("Familia", "Especie","Provincia","Localidad", "Fecha"),
rownames = FALSE,
options = list(
searchHighlight = TRUE,
language= list(url = '//cdn.datatables.net/plug-ins/1.10.11/i18n/Spanish.json'),
pageLength = 15))
primates <- inp[,c('species', 'eventDate')]
data <-primates %>%
group_by(especies = format(species)) %>%
summarize(suma_registros = n())
data$suma_registros<- as.numeric(data$suma_registros)
plot_ly(data, labels = ~especies, values = ~suma_registros,
type = 'pie') %>% config(locale= "es") %>%
layout(title = "Grafico Registros de Especies")
#Filtro de datos
AP<-filter(inp, species == "Alouatta palliata")
CC<-filter(inp, species == "Cebus capucinus")
AG<-filter(inp, species == "Ateles geoffroyi")
S<-filter(inp, species == "Saimiri oerstedii")
#Lectura de una archivo CSV con columnas de coordenadas
inp1 <-
st_read(
"primates-cr-registros.csv",
options = c(
"X_POSSIBLE_NAMES=decimalLongitude",
"Y_POSSIBLE_NAMES=decimalLatitude"),
quiet = TRUE)
# Asignación del sistema de coordenadas
st_crs(inp1) = 4326
# Mapeo
mapa <- inp1 %>%
dplyr::select(family,
species, stateProvince,
locality,
eventDate) %>%
leaflet() %>%
addProviderTiles(providers$OpenStreetMap.Mapnik, group = "OpenStreetMap") %>%
addProviderTiles(providers$Esri.WorldImagery, group = "Imágenes de ESRI") %>%
addCircleMarkers(
stroke = F,
radius = 2,
fillColor = 'red',
fillOpacity = 1,
popup = paste(
inp$family,
inp1$species,
inp1$stateProvince,
inp1$locality,
inp1$eventDate,
sep = '<br/>'),
group = "species") %>% addLayersControl(
baseGroups = c("OpenStreetMap", "Imágenes de ESRI"),
overlayGroups = c("primates")) %>%
addMiniMap(
tiles = providers$Stamen.OpenStreetMap.Mapnik,
position = "bottomleft",
toggleDisplay = TRUE)
#
alt <- getData(
"worldclim",
var = "alt",
res = .5,
lon = -84,
lat = 10)
altitud <- crop(alt, extent(-86, -82.3, 8, 11.3))
pal <- colorNumeric(
c("#006400", "#FFFF00", "#0000FF"),
values(altitud),
na.color = "transparent")
mapa %>%
addRasterImage(
altitud,
colors = pal,
opacity = 0.8) %>% addLayersControl(
baseGroups = c("OpenStreetMap","Imágenes de ESRI","Worldclim"))
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition